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The possibility to observe the neutrino-less double /3 de- 
cay and thus to prove the Majorana nature of neutrino as 
well as to provide a sensitive measure of its mass is a major 
challenge of to-day's neutrino physics. As an attractive alter- 
native we propose to study the inverse process, the radiative 
neutrino-less double electron capture Qv2EC. The associated 
monoenergetic photon provides a convenient experimental sig- 
nature. Other advantages include the favourable ratio of the 
Qv2EC to the competing 2v2EC capture rates and, very im- 
portantly, the existence of coincidence trigger to suppress the 
random background. These advantages partly offset the ex- 
pected longer lifetimes. Rates for the 0y2EC process are cal- 
culated. High Z atoms are strongly favoured. A resonance 
enhancement of the capture rates is expected to occur at en- 
ergy release comparable to the 2P—1S atomic level difference. 
The resonance conditions are likely to be met for decays to 
excited states in final nuclei. Candidates for such studies are 
4considered. The experimental feasibility is estimated and 
found encouraging. 

PACS numbers: 23.40.-s, 14.60.Pq, 23.40.Bw 

Introduction. The existence of massive neutrinos 
and the Dirac or Majorana nature of these particles are 
among the most intriguing topics of the present day 
physics. If neutrino is a Majorana particle then by def- 
inition it is identical to its charge conjugate. Thus the 
neutrino produced in one weak interaction vertex may 
be absorbed in another one. This leads to the nuclear 
reaction Qvf3(3: 

(A,Z-2)-> (A,Z) + e + e (1) 

(see fig.l). Amplitudes for such a process are propor- 
tional to the Majorana neutrino mass. While the ex- 
act value of the neutrino mass deduced from a successful 
0v(3(3 experiment is model dependent, the mere obser- 
vation of the effect proves unambiguously the Majorana 
nature of neutrino as well as the nonconservation of the 
lcpton number. This remains true regardless of the mech- 
anism causing the decay [1], [2]. 

The Oz//3/3 process proceeds via the emission of two cor- 
related (3 electrons. Its unique signature is that the sum 
of energies of the two electrons is equal to the total decay 
energy. The experiments searching for the §v(3~ f3~ de- 
cay can be divided into two categories: the calorimetric 
experiments, in which the material of the source is usu- 
ally identical with that of the detector, and the tracking 
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experiments, in which the source and the detector are 
separate. The former automatically sums up the ener- 
gies of the charged particles emitted. Large quantities 
of the material can be used. The main difficulty rests in 
suppressing the random background. The only way to 
attain this is by requiring extreme shielding conditions 
as well as the extreme purity of all the material of the 
detector housing and of the surrounding. The tracking 
detectors, counting the two electrons in coincidence, are 
somewhat less sensitive to the background. On the other 
hand, there are practical difficulties in handling the nec- 
cessarily large amount of the source material, of the order 
of tons, in the form of thin sheets sandwiched between 
the detectors. The detectors in both kinds of experiment 
must fulfill the high resolution requirement. Otherwise 
the Qvj3~ f3~ peak in the sum spectrum will not be dis- 
cernible from the dominating continuous physical back- 
ground due to the 2vf}~ (3~ decay. For the description of 
the problems involved we refer to reviews [2] - [6]. 

We suggest to study the inverse of neutrino-less double 
(3~ decay, i.e. the neutrino-less double electron capture. 
The excess energy is carried away by a photon. Crude 
estimates [7] for such a radiative process are encourag- 
ing, suggesting feasible experiments. These are discussed 
in the next section on the basis of a theory of the radia- 
tive capture adapted to the double electron capture case. 
The point of special interest is the resonant effect which 
happens when the photon energy equals the energy of 
atomic 2P — IS transition. In fortunate situations this 
effect may strongly enhance the OVcc rates. 

There are several experimental advantages of the ra- 
diative electron capture process : 

• the monoenergetic photon escapes easily from fairly 
thick layers of the source material without energy degra- 
dation; 

• the source can be separate from the detector; 

• the physical background due to the competing 2iseej 
process is low [8] ; 

• the photon emission is followed by that of the K X-ray; 
this provides a precious coincidence trigger to suppress 
the overwhelming random background; 

• low decay energies are favoured. Decays to excited 
states can thus be considered in realistic experiments in 
contrast to the (3f3 decays. The 7 transitions which follow 
offer yet another characteristic coincidence trigger. 

The price to pay is a sizeable reduction of the transition 
rates. The experimental questions and the optimal choice 
of the isotopes are discussed in the final sections. 

Amplitudes for the radiative capture. The 
dominant amplitudes for the radiative double electron 
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FIG. 1. Diagrams for the 0^ double beta decay and double 
electron capture processes. 



capture can be factorised into weak-nuclear, weak- 
leptonic and radiative factors [9], [11]. The expressions 
for radiative factors presented below are based on the 
theory for the radiative single electron capture process 
[12], [13]. 

The amplitude for Qv(3(3 decays from the {atom + 
nucleus) ground state | 0) to a final nuclear state Nf 
and final two electron state nf is 



00 _ 
Ov — 



(0\H w \n f ,N f ). 



(2) 



This differs from the amplitude for the double electron 
capture 



R^ = (n f ,N f \H w \0) 



(3) 



by the electron wave functions. Here fif denotes two 
electron holes in the final atom (say 15, 25 or 15, 2P). 
The reversed transition cannnot be tested in the same 
nucleus and Rq c v does not correspond to a real physical 
transition. Assuming the radiative process instead, we 
can write the corresponding amplitudes in the second 
order perturbation theory, to the leading order in the 



radiative interaction fl™, as 



FQ jv = E niiNi (nf,Nf | E 1 | ni,Ni) 



(m,Ni \H W \0) 
E - E(m) - E(Ni 



where H w is the weak interaction describing the neutrino- 
less process of fig.l. Intermediate states | ni,Ni) in the 
radiative process may involve ground and excited nuclear 
states. Only the electron is considered to radiate as the 
7 emission by the nucleus is unlikely [10], [9]. 

We look first for the transitions which conserve the 
nuclear angular momentum (0 + — > + ) as these offer the 
largest nuclear matrix elements. The final photon and 
its spin is generated by the electron radiation. Two basic 
processes are possible: 

(a) One of the electrons captured from the initial 
atomic state Ui„i radiates while it propagates towards 
the nucleus. 



(b) Two electrons are captured in a virtual process 
which generates a final atom with two electron holes. 
This final atom radiates and one of the holes is filled. 

Process (a) is standard in the single electron radia- 
tive captures. The description of electron propagation in- 
volves the electron Coulomb Green's functions [12], [13]. 
The two electron capture requires some modifications. 
To conserve the angular momentum, for one electron in 
IS state the other one must be in a higher nS state. The 
transition is of magnetic type. The electron wave func- 
tions must be antisymmetrized. 

Process (b) is less important in general. However, for 
small photon energy Q the virtual two electron hole state 
may be degenerate to the final atom-photon state (see 
below, eq. (12); see also refs. [15], [16], [17]). This leads 
to a singularity in Rq v . The most interesting situation 
occurs for Q = E(2P) — E(1S) i.e. when the final photon 
resonates with the 2P — 15 transitions in the final atom. 
The capture rate is enhanced up to the limit given by the 
natural K-L line widths. 

The cases (a) and (b) require different descriptions. 
In the simplest theory of Majorana neutrino [9], [11] the 
matrix element is: 

iC = 2 (J^j J d*dy{J N {*ym v h v {* - y)J N {y) a ) 

x*(n 1 ,x) 7p (l-75) 7cr * c (n 2 ,y). (5) 

The nuclear weak currents are denoted by J/v and the 
effect of neutrino propagation is included into the "neu- 
trino potential" with h v (r) w exp(~qr)/r, where q is 
an average momentum carried by the neutrino (the clo- 
sure approximation over nuclear states is used here). The 
projection on left-handed intermediate neutrino brings 
about the neutrino mass factor m„. The electronic part 
of this formula contains atomic wave functions \& and 
the charge conjugates \& . For the + — ► + transi- 
tions the nuclear part is reduced to Fermi and Gamow- 
Teller matrix elements defined by Mp = (0 | h v \ 0) and 
Mqt = (0 | h u a\U2 | 0). These enter via a combination 
M " = M GT - (^-) 2 M F and in this way the transition 



(4) matrix element (3) is brought to the form 



i^ = 2 fil) M 0v m v [*{ ni , 0)5(1, 2)*(na,0)U, (6) 

which resembles the standard (3(3 decay expression, [9]. 
The electron wave functions are required in the nuclear 
region. These are ^(n e , 0)u, where $(n e , 0) is the ra- 
dial part of the large component and u is the Dirac 
spinor. The spin matrix element is given by 5(1,2) = 
u(ni)(l + 75)^ (712). In eq.(6) the electron wave func- 
tion is to be antisymmetrised. For two IS electrons this 
involves antisymmetric spin zero combination of spinors 
which compensates an "inner antisymmetry" built into 
operator 5(1, 2) by the charge conjugation. 
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To describe the radiative capture of the type (a), one 
function W is to be replaced by some function "J^ which 
takes care of the photon emission and electron propaga- 
tion. In eq.(4), in coordinate representation, this func- 
tion is given by the expression 



At smaller Q the capture amplitude indicates an inter- 
esting structure. The electron propagator in eq.(7) has a 
pole at Q =| E(1S) - E(2P) |= Q res . To elucidate the 
atomic physics involved in this radiation process, ampli- 
tude (4) is presented in terms of the initial state 



, J-y(r, Tlini) — X ni 



(r | ni)(m | H 1 | n ini ) 
E-E(m)-E(0) 



(7) 



*i„i) =| 0) + S nii 



{ni,Nj \H W \0) 
E - E(m) - E(Ni) 



rH,Ni), (12) 



which involves the sum over continuum and discrete 
states of the electron, i.e. the Dirac Green's function 
in the external Coulomb field of the nucleus (Glauber 
and Martin [12], [13]). This solution is now implemented 
into eq.(4) to give the radiative amplitude 



E>7 

n 0v 



*(m,0)S(l,2,e,Q) 7 *(n2,0) 7 

(8) 



where the leptonic spins and the photon polarisation e 
enter the last term 



S(l,2,e,Q) 7 = 



2m e 



For the discussion of A(Q) and B(Q) we refer to ref. 
[13]. If Q w m e then A, B are close to unity for both the 
IS and 2S electrons. The spin conservation requires the 
two electrons to be in the symmetric spin and antysym- 
metric space combinations. This causes cancelations at 
small Z. Otherwise the effect of symmetrisation is small. 

The rate for radiative no-neutrino process is now 



27rdk 
(27r) 3 2fc 



S(k - Q) | 



(10) 



to be summed over photon polarisation and possible elec- 
tron pairs. The result is 

r°^= ^ 4 (M GT -(g) 2 M f ) 2 (m,/m e ) 2 
x | *(1S,0)*(2S,0) | 2 ^ pol (S(ni,n 2 ,e,Q) 7 ) 2 (11) 

where the summation over spin factor (S) is to be taken. 
This factor takes care of the angular momentum con- 
servation. The best possibility is the IS, 2S pair cap- 
ture accompanied by a magnetic photon transition. This 
rate is scaled by the S electron wave function at the nu- 
cleus. One has ^(nis, 0) 2 = (Zm e a) 3 f 2 /it where factor 
fg comes from the weak singularity of Dirac atomic wave 
function. Averaged over the nuclear volume it becomes 
f s = (2RZm e a) x - 1 with A = y/l - (Za) 2 and the nu- 
clear radius R. The Z & factor arises and thus one is 
interested in the heaviest possible atoms. There the rel- 
ativistic singularity enhances the rate further. 



where E — En(0) + E atom (0) is the initial energy of 
the system composed of the nuclear part Em and the 
atomic component E at0 m- Equation (12) describes the 
virtual transitions of the initial state | 0) under weak 
interactions which include virtual two electron capture 
states. Some of these states may be degenerate with the 
final atom + photon state to within the natural width of 
the K-L line. Rate enhancements of the order of 10 4 — 10 6 
may be obtained. In the single electron description by 
Glauber and Martin the electron propagator in eq.(7) has 
a pole at Q =\ E(1S, Z) - E{2P, Z) \= Q res {Z) where Z 
refers to the initial atom. In the many electron process 
considered here we use the energy levels for final(Z — 2) 
atom (see also ref. [17]). We obtain 



x(u( ni )[iA(Q)aexQ-B(Q)eQ}(l+ l5 )u c (n 2 )). (9) T°^(Q) 



r r (2P -> is) 



[Q-Qres(Z-2)} 2 + [T r /2f 



pec 1 2 



(13) 



where T r is the radiative width of the final two electron 
hole atom. We use the experimental energy values for the 
Z — 2 atoms corrected for the screening differences in the 
ionised atoms; the combined width is 3^(15) + T r (2P) 
similarly as for the hypcrsatcllite K h X-ray spectra. 

Experimental feasibility, optimal targets. 
The usual choice for targets in the double (3 decay is 
motivated by the phase space that favors large energy 
release ( roughly Q^p dependence). In contrast, the 
phase space dependence of the double electron capture 
eq.(ll) is rather weak and the radiative factors favor 
small Q. Two capture situations can be considered. At 
Q larger than the electron mass, the "magnetic" radia- 
tive capture of IS, 2S pair dominates. The calculated 
rates (Table 1) are based on values Mqt ~ 0.6/ fm and 
Mf/Mqt ~ —0.3 obtained in typical (3(30v reactions in 
a number of nuclei [18], [19]. These values give a crude 
estimate also for the Of cc processes, although some cal- 
culations indicate a larger coupling, [20] . 

At small Q the capture amplitude is given by the res- 
onant conditions. Several candidates for such a resonant 
capture may be found. Examples are given in Table II. 
The region of practical interest is limited to Q < 80 kcV. 
Such conditions are likely to be met in transitions to ex- 
cited states of final nuclei. The rate predictions are ham- 
pered by the poor knowledge of mass differences. Be- 
cause of the sharp atomic resonances a few keV mass 
uncertainty results in rates changing by several orders 
(cf. Tabic II). The atomic g.s. mass differences AM are 
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TABLE I. The radiative neutrino-less capture rates from 
the IS, 2S states [1/y], m v = 1 eV. 



Atom 



abundance % 



QjkeV] 



r(2S, is) 

~z — 'V> — 





15.84 


1628(5) 


2 




IfCd 


0.88 


241(7) 


1 


lO" 31 


IfSm 


3.07 


1730(4) 


1 


io- 31 


6 l 2 Er 


0.14 


1781(4) 


2 


10 -30 


7 fw 


0.12 


66(5) 


2 


IO" 31 


80 Hg 


0.15 


729(3) 


4 


10 -30 



taken from refs. [23], [24] and the Q values are AM re- 
duced by the two electron-hole excitation energies. The 
precision needed is well below 1 keV. 

The experimental signature for the double electron 
capture in the resonance conditions will be the KX - 
K h X-ray coincidence and, in the case of decays to ex- 
cited states, the triple coincidence with the gamma rays 
deexciting these states. The K h line has a measurable en- 
ergy shift. The experimental feasibility arguments have 
to include the decay rate and the cost estimates. Leav- 
ing the cost arguments aside and assuming 1 ton of the 
source material with 100 % isotopic purity and the cor- 
respondingly larger amount of a high resolution detector 
(e.g. a multilayered source-detector sandwich for X-ray 
detection combined with the medium resolution larger 7- 
ray detectors) it seems possible to design feasible exper- 
iments for the Of double electron capture process. The 
count rates expected depend strongly on the energy dif- 
ference with respect to the resonance value (see Table 
II). The 112 Sn isotope indicated in ref. [17] as the best 
choice, now with the the recent mass determination seems 
less profitable. Much higher decay rates can be expected 
for 152 Gd and 164 Er g.s. — > g.s. decays. However, there 
are no convenient experimental signatures in these tran- 
sitions. Best signatures are offered in the cases of decays 
to excited states (e.g. 112 Sn, 136 Ce or 162 Er) where there 
are characteristic high energy 7-rays in addition to the 
fc Q X-rays and the resonant transitions. 

The physical background. One advantage of the 
Of 7 process is a convenient ratio of the signal to the phys- 
ical background generated by the dominant v v 7 channel. 
Here, we briefly estimate this ratio defined as 



R 



S/B 



r°^(Q) 

T vu iN D 



(14) 



where Njj is the fraction of photons from the dominant 
ff7 decay mode emitted into the region from the end 
of the spectrum Q down to Q — D/2 and D is the pho- 
ton energy resolution. More details may be found in ref. 
[8] . For easier comparison the two-neutrino radiative rate 
is presented in terms of the no-neutrino radiative rate 
T 0l/J (k 1 ), given by eq.(ll) or eq.(13), as a function of 
the photon energy 



TABLE II. The resonant situations. The radiative neu- 
trino-less capture rates per year R/y, and per year and ton of 
the isotope R/y • ton, m v — 1 eV. The uncertainties of R/y • 
ton are due to lcr errors in the mass determination. Q - pho- 
ton energies, Q r - the resonant energies, Qx,Q-y - energies of 
accompanying X-ray, 7-ray photons, E* - excitations of final 
nuclei, [keV]. 



Atom 


112 qr, 
50 ™ 


136 p„ 
58 <^ e 


IfGd 


162 p 
68 t^r 


68 M 


18( 
74 


abnd.% 


1.01 


0.19 


0.20 


0.14 


1.56 


0. 


f.st. 


0j,lS,2P 


0+,lS,2P 


g.s.,lS,2P 


1 + ,1S,2P 


g.s.,2S,2P 


g-5 


int. st. 


0+ 1S.1S 


0+ 1S.1S 


g.s.,lS,2S 


1+,1S,1S 


g.s.,2S,2S 


g-s 


AM 


1919.5(4.6) 2418.9(13) 


54.6(3.5) 


1843.8(3.9) 23.3(3.9) 


14 


E* 


1870.9 


2315.4 




1745.5 






Q 


18.0(4.6) 


60(13) 


39.7(3.5) 


36.1(3.9) 


6(4) 


68 


Qr 


23.8 


33 


39.7 


46.6 


1 


56 


Q-Qr 


-5.8(4.6) 


27(13) 


0.04(3.5) 


-10.5(3.9) 


5(4) 


12 


Qx 


23 


32 


6.5 


46 


6.5 


55 


Q'f 


1253.4 


1496.9 




1665.1 








617.6 


818.5 




80.6 






R/y 


7. i(n 30 


2 • 10~ 30 


2 • 10~ 25 


1 • 10~ 28 


2 • 10~ 32 


3- 


R/y • t 


1 - 10~ 2 


10 _1 - 10~ 3 


10 3 - 10 _1 


1 - 10 _1 


10 _1 - 10~ 3 


3- 



J dLS(Q 



(15) 



where the phase space element dL 
Rn = (4TTM 2l/ )/(M°"m l/ ) is a dimensionless ratio of the 



and 



nuclear two-neutrino and no-neutrino matrix elements. 
Typical values M 2v w 1, M^ v w 1/fm follow from nu- 
clear model calculations, [3], [11]. To obtain R-s/b, ex- 
pression (15) is presented as an integral over the photon 
energy distribution T 1/1/7 = W(k^)dk 7 and the back- 
ground contribution becomes 



rQ 
y N D = / 

J Q—< 



Q-D/2 



WK(fc 7 )dfc 7 . 



(16) 



Formula W(k 7 ) = {Q - fc 7 ) c T , " /7 (k^)R 2 N /6(2rr) 4 follows 
from eq.(15). This indicates a cubic cut-off at the end of 
the spectrum due to the phase space. The ratio Rs/b 1S 
now given by the interplay of this cut-off and the energy 
dependence of r 0l/7 (A: 7 ). 

First, we consider the magnetic type transition related 
to the IS", 2S electron capture. For characteristic values 
Q = 1 MeV and m v = 1 eV one obtains T ^ /T vv i = 
5 • 10~ 5 . The photon energy resolution D — 3 KeV would 
yield very convenient Rs/b = 2-10 6 . The case of resonant 
electron captures is more involved due to the rapid energy 
dependence in r 0ly7 (fc 7 ) The result is plotted in Fig. 2. 
The Q < Q r es situation is very favorable and Rs/b is 
about 10 9 . For Q > Q res + T r the Rs/b f aiis down as 
[Q — Qres]~ 5 - At Q — Q res = 1 keV the ratio is still very 
convenient, Rs/b ~ 10 4 , but the conditions deteriorate 
quickly with the increasing energy separation. 

The random background (RB) varies with the source 
material. Crudely, the effcct/BR ratios are expected 
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FIG. 2. Schematic plot for the signal/physical background 
ratios: i?g/s(res)- continuous line for the resonant transition 
and R$/ Bijnagn)- dashed line for the magnetic transition. 
The rate dependence r°" 7 (res) vs (Q — Q res ) is shown in 
arbitrary units to indicate the resonance width. 

to be improved in comparison to the calorimctric 
/3~/?~experiments [5] by a factor between 10 and 1000, 
depending on the X-ray and the 7-ray coincidence trig- 
ger. 

We conclude that detecting the neutrino-less double 
beta process of any kind, be it (3~ , /3 + , or EC, will mean 
discovering the lepton number nonconservation and prov- 
ing the Majorana nature of neutrinos. It will also pro- 
vide a sensitive measure of the neutrino mass, thus giv- 
ing answers to some of the most urgent questions of to- 
day physics. This search calls for a massive, diversified 
experimental as well as theoretical effort. The radia- 
tive neutrino-less double electron capture, proposed in 
this work, is a valid alternative to double /3 ± emission, 
with several experimental advantages. The predicted res- 
onance rate enhancement makes this search plausible. 
Exact rate predictions require precise measurements of 
masses and likeways calculations of nuclear matrix ele- 
ments. The former are feasible with modern techniques 
[21], [24]. The latter are strongly assumption dependent. 
The consensus [22] is that they can be trusted to about a 
factor 3. One way to minimize this uncertainty is to carry 
out measurements for several nuclear species, preferably 
on both sides of the mass parabolae. Thus the double 
electron capture search can be considered not only as an 
attractive alternative to the double /3~. Assuming both 
kinds of experiment successful, it will also provide a much 
needed complementary piece of information. 
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